Exactly Solved Model for an Electronic Mach-Zehnder Interferometer 
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We study nonequilibrium properties of an electronic Mach-Zehnder interferometer built from integer quantum 
Hall edge states at filling fraction iz—l. For a model in which electrons interact only when they are inside the 
interferometer, we calculate exactly the visibility and phase of Aharonov-Bohm fringes at finite source-drain 
bias. When interactions are strong, we show that a lobe structure develops in visibility as a function of bias, 
while the phase of fringes is independent of bias, except near zeros of visibility. Both features match the results 
of recent experiments [Neder et al. Phys. Rev. Lett. 96, 016804 (2006)]. 

PACS numbers: 71.10.Pm, 73.23.-b, 73.43.-f, 42.25.Hz 



Questions about phase coherence in interacting quantum 
systems out of equilibrium are of fundamental and wide- 
ranging importance. Despite great progress over the past 
decade, many aspects of nonequilibrium problems remain 
poorly understood. One recent example of this situation is the 
"unexpected behaviour" observed in state-of-the-art experi- 
ments on electronic Mach-Zehnder interferometers (MZIs) 
(d HI [stl driven out of equilibrium by an applied bias volt- 
age. In these experiments the visibility of Aharonov-Bohm 
(AB) fringes in the conductance shows a lobe-like structure 
as a function of bias, while the phase of oscillations is inde- 
pendent of bias even with different interferometer arm lengths, 
except at zeros of the visibility where it jumps by vr. 

These observations have attracted a lot of attention. It was 
immediately appreciated [1] that they lie outside a single- 
particle description. Moreover, since integer quantum Hall 
edge states scale to non-interacting chiral Fermi gases at 
low energy, the finite-range of electron-electron interactions 
seems to be crucial. The effort to understand interaction ef- 
fects in MZIs at integer filling is therefore linked with work on 
non-linear effects in non-chiral Luttinger liquids [4], as well 
as to interferometry of fractional quantum Hall quasiparticles 
Jstl- The most obvious consequence anticipated from interac- 
tions is debasing. This may arise from external noise |@] or 
internally [7, 8], but in both cases is expected to suppress AB 
fringe visibility smoothly with increasing bias, in contrast to 
observations. It has been found, however, that zeros in visibil- 
ity can arise if the edge channels that form the interferometer 
arms are coupled to another channel: such an extra channel 
may be a feature of sample design and is present intrin- 
sically at 1^=2 [To*]. Although those results are encouraging, 
they do not seem sufficiently universal to explain all current 



experiments. In this context, two recent papers 111 IL 11211 that 



obtain visibility oscillations from calculations of interaction 
effects at v—l represent an interesting advance. These papers 
contain illuminating physical insights, and similar phenom- 
ena have been shown to exist in another context (l^, but ap- 
proximations used in fill |T2ll are not standard ones and their 
reliability is hard to judge. 

In this Letter we present an exact calculation for a simpli- 
fied model of an interferometer. It reproduces the main signa- 



tures observed experimentally and shows that the lobe 

pattern is a many-body effect, which would not appear in any 
approximation that treats single particles moving in a static 
mean-field potential. The model is illustrated in the inset to 
Fig. [T] As in previous studies, two quantum Hall edge chan- 
nels, both with the same propagation direction, are coupled 
at two quantum point contacts (QPCs). The simplifying fea- 
ture of the model is that electrons interact only when they are 
inside the interferometer. This allows us to combine a descrip- 
tion of the contacts using fermion operators with a treatment 
of interactions using bosonization. Within the MZI we take 
interactions only between two electrons on the same arm and 
with fixed strength independent of distance, although it would 
be feasible to relax these restrictions. We consider an initial 
state in which Fermi seas in the two channels are filled to dif- 
ferent chemical potentials, to represent the bias voltage, and 
evolve this state forward in time using the Schodinger equa- 
tion. At long times the system reaches a stationary regime. In 
this regime we calculate current and differential conductance 
as a function of chemical potential difference and enclosed 
AB flux. Our main results are presented in Figs.[T]and|2] and 
discussed following an outline of their derivation; details will 
be presented elsewhere 01411 . 

The solution we describe is significant more broadly as a 
rare example of a solved non-equilibrium scattering problem. 
One earlier instance is that of tunneling between fractional 
quantum Hall edge states lITsll . while another is the interact- 
ing resonant level model, treated recently by a form of Bethe 
Ansatz [16], and using boundary field theor vjl 1711 . The re- 
markable structure observed experimentally UrUlsl] makes 
the MZI particularly interesting in this context. 

The Hamiltonian H — Hkin + Hint + Htun for the model 
has three contributions, representing respectively: kinetic en- 
ergy, interactions, and tunneling at contacts. We formulate H 
initially for edges of length L with periodic boundary condi- 
tions, then take the limit L ^ 00. Then 
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ri=l,2 



■ip'l{x)dx'>Pn{x)dx, (1) 



where vp is the Fermi-velocity and rj — 1, 2 is the channel 
index. The Fermi field operators can be written as (x) — 



2 



L J2k ^kriG^'^^ , with k = 2'iTnk/L and Uk integer, and 
{cfcjj, c^,,'} — 5kq5riri' ■ Interactions are described by 



Hi, 
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-L/2 



Uri{x,x')pn{x)pri{x')dxdx' , (2) 



where (x) = {x)i)ri{x) is the electron density operator 
In our model [/^(x, x') = for x, x' ^ (0, d,,). Finally, the 
QPCs are represented by 

^tu„ = fae*"Vi^(0)V'2(0) + V6e^'3^+(di)V'2(d2) + h.c. (3) 

The AB-phase appears here as ^pab = (i ~ ol. 

The total current / from channel 1 to 2 has contributions la 
and /b arising from each QPC, which can be written in terms 
of expectation values of operators acting at points infinitesi- 
mally before the QPC. Each contribution can be separated into 
a term that is not sensitive to coherence between the edges, 
and another that is sensitive. We define taf, — smQaft and 
Taf) — cos6'a,6 with da.h = Vafilfivp, and denote expectation 
values by (. . .). A straightforward calculation yields for QPC 
6 the expressions Ih = l\P + Z^'^'', with 



= evFtl{h{d^) ~ P2[d2)) 



(2) 
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h.c.l 



where G12 — V'j''(rfi)V'2(rf2)- Terms in are obtained from 
these for Ii, by replacing di and d2 with 0, and vi, with Va- 
Since there is no coherence between channels before QPC 

a, la =0 and the term responsible for AB oscillations in 

(2) 

current is . The bias voltage is = {fii — P2)/e and 
the differential conductance is C/ = edi/dpi (with p2 fixed). 
Q oscillates with lpab, having maximum and minimum val- 
ues Gma.x and Gmin, and AB fringe visibility is defined as 

(^max ^min)/ (^max ^min)- 

The central task is therefore to calculate the correlator 
(G12), and our approach is as follows, (i) We work in 
the interaction representation, evolving operators with Hq = 
Hkin + Hint and treating Htun as the 'interaction'. Then 
-0^ (x, t) = e'^o*/'^Vi^ (x) e-*^o*/'* (note that we distinguish 
operators in the Schrodinger and interaction representations 
by the absence or presence of a time argument). The wave- 
function of the system, denoted at t—Q by \Fs), evolves with 
the S-matrix S {t) = Texp{-(i/n) Htun (f) dt'}, where 
T indicates time ordering, (ii) Time evolution of operators is 
calculated using bosonization to diagonalise Hq. (iii) Results 
are written in terms of operators in the Schrodinger picture, 
with boson operators re-expressed using fermion ones. This 
yields an expression for G12 suitable for straightforward nu- 
merical evaluation. We next outline these three steps. 

Step (i): Evaluation of S{t) hinges on our restriction of 
interactions to the interior of the MZI. Specifically, sepa- 
rating Htun into parts H^^n and Hf^^n due to each QPC, 
we find from step (ii) that /ff\,„(t2)] = and 

[Gi2{ti), H'^{t2)] — 0, provided ti > t2- The first commuta- 
tor allows us to factorise the S-matrix as S{t) — S''{t)S"'{t), 



where S'^{t) is calculated using -fff„„ and S^{t) using Hf^n- 
The second ensures that [S''^{t)]+Gi2it)S''{t) ^ Gi2(t), so 
that an explicit form for S^{t) is not required in the cal- 
culation. Since QCP a acts before interactions, S°'{t) is 
easy to evaluate: we have 7?"jj„(t2)] = for any 

ti,t2 > and so may omit time ordering. In particular, we 
will need to compute the action of S'"(t) on fermionic opera- 
tors. It is a rotation in the space of channels and can be written 

^Pnix) = [S''{t)]+i;rj'{x)S''{t). ForO < X < vptwe find 
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Step (ii): We compute time evolution under ifo using 
bosonization jlSll . Fermion operators are written in the form 



(5) 



where Fj, are Klein factors with commutation relations 
{Fr,^ F^,} ^ 2Srin' ^'^'^ bosonic fields are defined as 

,-91/2 

q>Q 



(^) = - E (W'/i)'^' {e'''%rj + h.c.)e-w^^ (6) 



with a an infinitesimal regulator Plasmon creation operators 
obey bosonic commutation relations [bq,-; , b'^,y] = SqkS,pf and 
are expressed for q > in terms of fermions as 



k— — oo 



(7) 



Since Hq does not couple channels, we restrict attention 
to a single channel and omit channel labels until we reach 
step (iii). The kinetic energy Hkin for a single edge has the 
bosonized form 



Hkin — 



dx 



2 .'-L/2 27r 



1( ^^2 , 'i.TThvF 

{da,(t)ix)) — 



N{N + 1) (8) 



where N = c^Ck is the particle number operator Simi- 
larly, Hint is quadratic when written using the bosonic repre- 
sentation of the density operators, p (x) = ~^dx(l){x) + 
N/L. The time dependence of (f>(x,t) can be found by 
solving the equation of motion. Since our choice of non- 
uniform interactions leads to a coupling between the plas- 
mon and number operators, we make the separation </> (x, t) = 
(^(°)(x,i) + <^(i)(x,t), where 0(o)(x,i) oc N/L and 
(f)^^'' (x, t) is independent of N, satisfying 

2Trh{dt + VFd,)4>^^^ (x, t)^-J U{x, y)dy4>^^^ {y, t) dy . 

(9) 

The solution can be written in the form 

0(1) (x, t) = K{x, y- 1) [4> (y) - 0(") {y)]dy , 

J -L/2 
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FIG. 1: Visibility as a function of bias voltage for MZI with di — d2 
and to = t\ = 1/2 at interaction strengths: 7 = 1 (dot-dashed 
line), 7 = 2 (dashed line), 7 = 3 (full line), where 7 = 27r7. Inset: 
schematic view of model studied. 

where the Green function K{x, y; t) can be constructed in the 
usual way from the eigenfunctions of the time-independent 
equation, 

27rhvFidx -ip)fp{x) ^ - J U{x,y)dyfp{y)dy. 

We now specialise to interactions that are constant within 
the interferometer: U{x,x') — g for x,x' G (0,0?) and 
U{x,x') = otherwise. This form of the potential is the 
one treated approximately in fll']. It is characterised by the 
dimensionless coupling constant 7 = gd/2TThvF- We find in 
the limit L ^ 00 

( e*P^ a; < 

fp (x) ^ } rp + Spe^'P'-^ < X < d . 
[ e'P^'-'^p x>d 

The coefficients Sp — {1 + tp)^^ and Vp — tpSp, with tp = 
{ij/pd){l — e'P''), are obtained from matching fp{x) at a; = 
0, d. The phase shifts of plasmons 6p due to the interactions 
are given by e~"'p = (1 + t*)/{l + tp). Similarly, we find 

^(0) (x) = 2'K^Nx/Lforx G (0,d), where 7 = 7(1 + 7)^1. 

In this way we find an expression for K{x,y;t). Setting 
X = d,it simplifies at long times to 

1 f°° 

K{d,y;t) = — / dpe'^P^'^-y-"^*^-^'-^ . (10) 

Using this and Eqns. ^ and (|2l), we write 0*^^^ (x, t) as a bi- 
linear in the fermion operators and Ck ■ 

Step (iii): We employ this result to construct an expression 
for (?i2 {t) in terms of fermion operators in the Schrodinger 
representation. To this end, we start from Eq. (|5]l in the in- 
teraction representation at time t and substitute for 0,, (d,, , t) 
as described. We also eliminate the combination T = 
(27ra)^^/^F^e*T^'''^'' by inverting the bosonization identity, 
Eq. (|5]l, writing 



for = dfj — vpt. Finally, we substitute for and 
in 4)ri{x) using Eq. (|7]i. The result (omitting an unimportant, 
constant phase) is the operator identity 

i;^{d^,t)=e-"^-^^{z^) . (11) 

Here Q,, = I^oc^^^''^ ~ z,^)p,^{x)dx, where the kernel 

Q,-i{x) = L^^ X]^-oo Qnil)^^'^^ ^'^^ f^i" choice of in- 
teraction the Fourier transform 

Q,M) - 27r^d,foiQdr,/2){l+^e~''''-/'Mqd,/2))-\ 

(12) 

in which jo (a;) — x^^ sinx. 

In this way we arrive at the expression 

{G,2{t)) ^ e''^{Fs\[S-{t)]+4>t{zi)e'H2iz2)S-it)\Fs) . 

Here $ is an initial phase that is independent of voltage, and 
R ^ Qi- Q2- The action of S"^(<)+ and S'^it) on the op- 
erators they enclose is given by Eq. (|4|l, and evaluation of 
{Gi2{t)) reduces to the calculation of correlators of the form 
C^r, = (Fs|c+exp(iX;Q/3Ma^c+C0)c^|Fs), where the in- 
dices specify both channel and momentum, and the matrix 
M is obtained from [S°- {t)]+ RS'' {t) . One can show that 
C^n = D,^^ det D with D constructed from the matrix ele- 
ments of exp(iM) between the single-particle states that are 
occupied in the Slater determinant \Fs). We calculate C^^ 
numerically, achieving convergence of the results when keep- 
ing up to lO'' basis states and 400 particles in each channel. 

The physical interpretation of the solution we have pre- 
sented is as follows. Each electron passing QPC b at time 
t has an accumulated phase from its interactions with other 
electrons. The phase is a collective effect and is represented 
by the operator Q,, in Eq. ( fTTT) . Contributions from interac- 
tions with particles at a distance x from the one at QPC b 
have a weight determined by the kernel Q,-i{x), illustrated in 
the inset to Fig. |2] This weight is largest near x = 0, show- 
ing that interactions with nearby electrons are most important. 
Moreover, since — for x < — d^, a given electron is 

uninfluenced by the ones behind, that enter the interferometer 
after it exits. The precise form of the kernel reflects the full 
many-body physics of the problem: a similar kernel appears 
in Eq. (11) of Ref. Il ill , but with a simpler form because of 
the approximations employed there. 

A consequence of the phase Qn is that many-particle in- 
terference influences the MZI conductance. As an illustra- 
tion, consider the quantum amplitudes for two particles to 
pass through the interferometer on all possible paths connect- 
ing given initial and final states. Paths for which both particles 
propagate on the same arm of the interferometer have an inter- 
action contribution to their phase that varies with their sepa- 
ration and is absent if the two particles propagate on different 
arms. Destructive interference between paths with different 
interaction phases generates the observed lobe structure. 

We now turn to our results. The parameters in the model 
are: the dimensionless interaction strength 7, the transmission 
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FIG. 2: Dependence of interference fringe phase on bias voltage 
for MZI with unequal arm lengths, d2/d\ = 1.2, at interaction 
strengths: 7 = (dot-dashed line), 7 = 3 (dashed line), and 7 = 10 
(full line). Inset: the kernel Q{x) of Eq. l[TT]l at 7 = 10 (full line) 
and that of Ref. 1 1 1] (dashed line). 

probabilities and tl, the ratio d2 / di of arm lengths, and the 
dimensionless bias voltage eV^/did^/iTrhvp. We consider 
1 < 2ttj < 10, 1 < c?2/c?i < 1.2 and first discuss behaviour 
with = tl = 1/2. The dependence of visibility of AB 
fringes on bias voltage and interaction strength is presented in 
Fig. \T\ taking equal arm lengths and transmission probabili- 
ties of 1/2 at both QPCs. The key features of all three curves 
in this figure match those of the experiment (see Figs. 2 and 
3 of fT]): with increasing bias there is a sequence of lobes in 
the visibility, which have decreasing amplitude and are sepa- 
rated by zeros. The phase of AB fringes is also influenced by 
interactions. Results are displayed in Fig|2] For an MZI with 
different arm lengths (as in this figure), the fringe phase with- 
out interactions varies linearly with bias, because the Fermi 
wavevector fc^ is linear in bias and the phase difference be- 
tween particles traversing the two arms is kF{d2 — di). With 
increasing interaction strength the phase dependence on bias 
develops into a series of smooth steps, each of height tt. The 
risers of these steps coincide with minima of the visibility. 
Strikingly, with strong interactions phase steps at minima of 
the visibility persist for di = d2, even though in this case 
phase would be independent of bias without interactions. The 
stepwise phase variation we find at large interaction strength 
also matches observations (see Fig. 2 of 1 1]). 

Behaviour is insensitive to the transmission probability at 
QPC b, apart from the overall scale for visibility. Departures 
from t'^ = 1 /2, however, eliminate the exact zeros in visibil- 
ity, leaving only sharp minima. A difference in arm lengths 
has a similar though much weaker effect. 

The width in bias voltage of the central visibility lobe de- 
fines an energy scale. In our model this scale is of order g at 
large 7. Taking vp = 2.5 x lO^'ms^^, d ~ 10/im and the 
permittivity e ~ 12.5 of GaAs, we estimate from the capac- 
itance of an edge channel g ~ 10/ieV. This is similar to the 
experimentally observed value of about 14 fieV 



Our calculations rely on a simplified form for interactions, 
but we believe our choice is quite reasonable. Our central 
approximation is to neglect interactions between an electron 
inside the MZI and one outside. In practice, such interactions 
will anyway be screened by the metal gates that define the 
QPCs. We also neglect interactions between a pair of elec- 
trons that are both outside the MZI. This is unimportant: be- 
fore electrons reach the MZI, such interactions do not cause 
scattering because of Pauli blocking, while after electrons pass 
through the MZI, these interactions cannot affect the current. 
Within the MZI we represent interactions by a capacitative 
charging energy. Such a choice is standard in the theory of 
quantum dots and has been applied previously to interferom- 
eters 13 El- 

In summary, we have calculated the visibility of Aharonov- 
Bohm fringes in the differential conductance of an electronic 
MZI out of equilibrium, taking exact account of interactions 
between electrons. From our calculations we obtain a lobe 
pattern in the dependence of visibility on bias, and jumps in 
the phase of fringes at zeros of the visibility, as observed ex- 
perimentally [1, 2, 3]. 
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